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GLLH EM Invisible Cloak With Novel Front Branching And 
Without Exceed Light Speed Violation 

Ganquan XicI and Jianhua Li, Lee Xie, Feng XiqJ 

GL Geophysical Laboratory, USA 

(Dated: July 1, 2011) 

In this paper, we propose new Global and Local (GL) electromagnetic (EM) cloaks with distinctive 
class material Qq/j [og°' (ba$/h)h^ (GLLH Cloak) without exceed light speed violation. The refractive 
index of the GLLH cloak material, n{r), is large than one or equal to one. Our GLLH cloak is 
created by GL EM modeling and GL EM cloak inversion with searching class a^ 13 log" {ha 13 / h) h^ . 
The GLLH cloaks in this paper have finite speed and have no exceed light speed? physical violations 
and have more advantages. The GLLH EM cloaks can be practicable by using normal materials 
and are available for all broad frequency band. The GL EM cloak inversion and electromagnetic 
integral equation for cloak are presented in this paper. The novel EM wave propagation and front 
branching in the GLLH cloak by GL EM modeling are presented in this paper. The EM wave front 
propagation in GLLH cloak is behind of the front in free space. At time steps 118di, in the GLLH 
cloak, the wave front is curved as a crescent like and propagates slower than the light in free space. 
At the time step 119df , the EM wave inside of the GLLH cloak propagates slower than light speed, 
moreover, its two crescent front peaks intersect at a front branching point. At the front branching 
point, the front is split to two fronts. The novel front branching and crescent like wave propagation 
are displayed in the following figure 1, figure 2 and figures 5 -20 in this paper. All copyright and 
patent of the GLLH EM cloaks and GL modeling and inversion methods are reserved by authors in 
GL Geophysical Laboratory. 

PACS numbers: 13.40.-f, 41.20.-q, 41.20.jb,42.25.Bs 



I. INTRODUCTION 

In this paper, we propose new Global and Local (GL) 
electromagnetic (EM) cloaks with distinctive class mate- 
rial ttajs log" {bai3/h)h^ (GLLH Gloak) with finite speed, 
and without exceed light speed violation. The refrac- 
tive index of the GLLH cloak material, n{r), is larger 
than one or equal to one. Where Uap is cocfhcients, 
h = r — Ri. Our GLLH cloak is created by GL EM mod- 
eling and GL EM cloak inversion with searching class 
aap log" {ba 13 / h)h^ . We name our cloak as GLLH cloak. 
The GLLH EM cloaks can be practicable by using nor- 
mal materials and are available for all broad frequency 
band. In 2001, we used Global Integral and Local Differ- 
ential (GILD) EM modeling and inversion to detect the 
fly model imaging from EM sources in some frequency 
band. A strange double layer clothes phenomenon was 
discorvered in the residual field. The double layer cloth 
protects the fly object from the exterior EM fleld de- 
tection. We called the strange phenomena as GILD ef- 
fect and published it in paper [1]. In this paper. We 
present the double layer clothes imaging again in Figure 
3 and Figure 4. For investigating the strange dou- 

ble layer clothes phenomena, we developed Global and 
Local fleld electromagnetic modeling method [2] and GL 
Metro Carlo inversion method [3]. These new methods 
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are complete different from the traditional Finite differ- 
ence scheme [4] and Finite Element Method [5] . The GL 
method has totally advantages than the traditional meth- 
ods. The big matrix solving cost and error reflection on 
the artificial boundary are difficulties in the traditional 
numerical simulation methods. The GILD EM model- 
ing and inversion [6] are important progress to overcome 
these difficuhies. The GL EM modehng and GL Metro 
Garlo inversion have advantages to completely overcome 
these difficulties in history [2] and [3]. In the physical 
and quantum mechanical simulations. Born approxima- 
tion [7] is often used. However, the Born approxima- 
tion is not accurate and can not be used in the high 
frequency band and high contrast and singularity mate- 
rials. We used the GL EM modeling and inversion to 
simulate the double layer clothes phenomena [1] and dis- 
cover a mirage image [8], and proposed GL double layer 
cloak [9-13]. Pendry at al proposed the single layer cloak 
by transform optical technology in [14] (Ps cloak). How- 
ever, the exceed light speed and infinite EM field speed 
are its two main physical violations. In [15], we proved 
a theorem that there exists no Maxwell wavefield can be 
excited by sources inside the single layer cloaked con- 
cealment, if the cloak material parameters are finite, and 
the concealment is fill by basic EM materials eo and /io 
. By the theorem, there is no EM wave can be excited 
in Ps' cloaked concealment that is its thrid physical vi- 
olation. Our GL double layer cloak [9-13] overcame this 
difficulty. The GL outer layer cloak has invisibility func- 
tions. The GL inner layer cloak decays the internal wave 
field, such that the internal wave field can not propagate 



to outside of the inner layer. There is no exceed light 
speed violation in the GL inner layer cloak. In the GL 
outer layer cloak, we obtained reduced degenerate cloak 
material, when r — >■ i?i er^J■e ~ {r — Ri) — >■ 0. How- 
ever, Ps cloak is strong degenerate material, in which[14] 
the Erfie = R^r - Rif/r'^/{R2 ~ Rif , when r ^ Ri 
ErUs « (r — i?2)^ — > 0.. We proposed a GLWF dou- 
ble layer EM cloaks in broad frequency band [16] to 
overcome the exceed light speed. The Figure 9 in [16] 
shows that the EM wave propagation in Ps cloak does 
exceed light speed . The Figure 7 in [16] exhibits that 
the EM wave propagation in GLWF double layer cloak 
does not exceed light speed. As a great jump of the 
GLWF cloak [16], by using GL Metro Carlo inversion 
[3] , we propose the new GLLH cloaks with class material 
Oq/S '^og" {ba is/ h)h^ , one of which is formulized by (2)-(5) 
in this paper. The cloak has nonzero and positive e^, fj-r, 
eg, and /ig in whole cloak domain Ri < r < Ro, that is its 
large advantage over than the 2rd order strong singular 
Ps cloak with zero parameters on the inner boundary, 
Er = and Hr = 0. We proved that the refractive index 
of our GLLH cloak is larger than one or equal to one, 
n{r) = y'Brl-J'e > 1, in whole cloak domain Ri < r < Ro, 
when R, = 0.249, Ro = 0.47. By using GL EM modeling 
[2] , the novel full EM wave propagation in the cloak show 
that there is no exceed light speed violation in our GLLH 
cloak in in this paper. The novel EM wave propagation 
and distinct front branching in the GLLH cloak by GL 
EM modeling are presented in in this paper. The EM 
wave front in GLLH cloak is behind of the front in free 
space. The figures 8-13, show that the wave front is suc- 
cessively curved as a crescent like and propagates slower 
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FIG. 2: (color online) Electric wave propagation at time step 
119dt, the two peaks of the crescent like wave front intersects 
at a branching point, red S in right is source location 



than the light in free space. In particular, at time steps 
ll8dt, in the GLLH cloak, the figure 1 shows the curved 
EM wave front is much similar with the crescent. At the 
time step 119dt, the electric wave inside of the GLLH 
cloak propagates slower than light speed, moreover, its 
two crescent peaks intersect at a novel front branching 
point. At the front branching point, the front is split 
to two fronts. One is outgoing front, which is propagat- 
ing forward to out of the cloak and bring out most wave 
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FIG. 1: (color online) Electric wave propagation at time step 
118dt, the wave front is curved as the crescent like, red S in 
right is source location 



FIG. 3: (color online) A novel double layer cloak surround 
space shuttle model to prevent detection has been discovered 
in GILD EM inversion in 2001, which was published in SEG 
expand abstract in 2002 [] . 



energy. After the outgoing front propagating out of the 
GLLH cloak, it is recovered to the original wave front in 
free space. Other front is similar with a sphere surface 
wave and is attractive propagating to the inner bound- 
ary. The attracting front propagation is very shower. Its 
amplitude and speed are rapidly decay to zero. The at- 
tracting front is propagating to inner boundary, but it 
can not reach to the inner boundary r — Ri. The novel 
front branching and crescent like wave propagation are 
displayed in the figure 1 figure 2 and figure 5 -20 in this 
paper. 

When source is located outside of the GLLH cloak, the 
observer point r is located inside of the cloak and it is 
going to the inner boundary r — Ri, the EM wave field 
delay to zero inverse radial. On the inner boundary, the 
radial electric wave field Er{Ri, 0,4)) =0 and the radial 
magnetic wave field Hr{Ri, 0, 0) = are distictive prop- 
erty of the GLLH invisible cloak. On the inner boundary 
of other transformed cloak, the radial wave EM field are 
not zero. There is a novel geometry for our GLLH invis- 
ible cloak. 

Ulf Leonhardt et al. proposed a new cloak with fi- 
nite speed based on a Euclid and non Euclid joint trans- 
form[17]. The ULF cloak overcomes the "infinite speed" 
physical violation, even though its refractive index less 
than one in some subdomain, for example, n((T, cr',r) = 
ri(0.757r,7r, r) < 1. Chen et al and Zhang et al proposed 
an analysis for Ps cloak using Mie transform [24] [25] . 

The creating motivation of our GLLH clock is accord- 
ing to our double layer cloth observation in 2001 which is 
shown in Figure 3 and Figure 4, and our 3D FEM super 
convergence observation in 1973. [5] [18] [19]. 

The form Uap \og'^ {hap / h)h^ is used as estimation form 
in number theory [20] [21] [22], PDE, numerical mathe- 
matics [5] [23], probability, geometry, and in other estima- 
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FIG. 4: (color online) A novel double layer cloak surround 
fly model anti detection has been discovered in GILD EM in- 
version in 2001, which was published in SEG expand abstract 
in 2002 [1 . 



tion theory in physics and other science. However, no 
one think it can be used in invisible cloak material be- 
fore this paper. Fortunately, the form aap \og"^ {hap / h)h^ 
is used as search class of the GL cloak inversion to cre- 
ate the GLLH invisible cloak that is our inspiration and 
discovering. The finding inversion is to find object from 
the observer measured scattering wave from the object 
and is excited by control or natural sources. The invis- 
ible cloak generation is hiding inversion. The invisible 
and hiding inversion is to generate a cloak material such 
that the scattering wave from the cloak is zero. The dou- 
ble layer cloth and fiy imager in Figure 3 and Figure 4 
show that the finding inversion and hiding inversion are 
contradiction and compromise each other. 

The description order of this paper is as follows. In- 
troduction is described in Section 1. In Section 2, we 
propose new GLLH em cloak without exceed light speed. 
The GL EM cloak inversion is proposed in Section 3. 
In Section 4, we present novel wave front branching and 
EM wave propagation in the GLLH cloak without exceed 
light speed. The distinct invisible properties of GLLH 
EM cloak are proposed in Section 5. In Section 6, we 
propose GLLH EM double layer cloak. In section 7, we 
present that the EM wave field decay to zero inverse ra- 
dial. In section 8, we describe history and discussion. 
The conclusion of this paper is described in section 9 



II. GLLH EM CLOAK WITHOUT EXCEED 
LIGHT SPEED 

In this section, we propose a new GLLH EM invisible 
cloak for broad frequency band and without exceed light 
speed violation. 



GLLH EM INVISIBLE CLOAK 



In the Maxwell equation 



V xE^-fl§H + M 

V X iJ = e§E + J. 



(1) 



Where E is eletric wave vector field, .ff is magbetic wave 
vector field 



£ = diag [ Sr £e £<t> ]eo,Ri <r < Ro, (2) 

p. = diag [ fir f^e A^^ ] Mo, -R* <r<Ro- (3) 

We propose the GLLH EM invisible cloak material pa- 
rameters in 3D spherical coordinate system, the cloak 
material parameters are radial dependent as follows. 



2Ri{Ro-Ri 



^^^ ^ l+4Ra+2Rl y 
2+6Ra 
"log(ei/-R.o(_Ro__R„)/(r-fli)) 

1+2R„ (r-fii) 



I 2+6Ro (r--Ri) 



(4) 



ee = e^ = /ig 



1 



1 
"21og3(ei/Ro(_Ro-_R,)/(r-_R,))(r-fli) 



(5) 



where the Ri < r < Ro is the spherical annular cloak 
domain, Ri is inner radius of the annular cloak, Ro is 
outer radius of the annular cloak, e is the dielectric pa- 
rameter matrix, p, is the permeability parameter matrix. 
So is the basic dielectric parameter, /xo is the basic mag- 
netic permeability parameter, e^ is relative dielectric in r 
direction, ee is relative dielectric in 6 direction, e^ is rel- 
ative dielectric in cj) direction, fir is relative permeability 
in r direction, fig is relative permeability in 9 direction, 
ficjy is relative permeability in (j) direction. 



A. 3D EM INTEGRAL EQUATION 

A 3D electromagnetic integral equation is proposed in 
this section as follow: 



E{r,r,,t) 



Eb{r,rs,t) 
Hb{r,rs,t) 



(6) 



2. function II 

When the source r^ is located outside of the cloak and 
observer r is inside of the cloak , we verify the invisible 
function II of the GLLH EM cloak materials in (4)-(5) as 
follows: 

Exterior EM field does not penetrate into concealment 
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B. INVISIBLE FUNCTIONS OF GLLH EM 
CLOAK 

1. function I 

When the source r^ and observer r are located outside 
of the cloak, we verify the invisible function I of the 
GLLH EM cloak materials in (4)-(5) as follows: 

Exterior EM field must not be interfered by scattering 
from the cloak : 



(9) 

Where E{r,rs,t) is the electric intensity field, 
H{r,rs,t) the is magnetic intensity field, r, is observe 
space location variable, r^ is source space location vari- 
able, t is time, Eb{r,rs,t) is the background electric in- 
tensity field, Hb{r,rs,t) is the background magnetic in- 
tensity field, G^^^(r?r, i) is 6 x 6 EM Greens tensor 
matrix, diag onal [5 Dn, 5 D 22] is is 6 x 6 variation EM 
material matrix as follows: 



1 = (^- '^0l)gt, 01^22 



SDi 

e = diag [ £r 
j2 = diag [ fi 



d 



,SD2 

£e £4, ] eo, 
Me A^0 ] A^o, 



(m 



Mo/)^ 



dt 



(10) 



Where e is the dielectric parameter matrix, p, is the per- 
meability parameter matrix, the dielectric and perme- 
ability can be isotropic or anisotropic materials, eo is the 
basic dielectric parameter, ^0 is the basic magnetic per- 
meability parameter, Sr is relative dielectric in r direc- 
tion, ee is relative dielectric in 9 direction, e^ is relative 
dielectric in cf) direction, fir is relative permeability in r 
direction, fig is relative permeability in 9 direction, fi^ is 
relative permeability in cf) direction, for the cloak mate- 
rial, the e and fi are proposed in formulas (3) — (5) in 
this paper. 



E{r,rs,t) 
H{r,r,,t) 



= 0, 



(7) 



Our GLLH EM cloak completely overcomes the exceed 
light speed?and infinite speed? two physical violations. 
The GLLH EM cloak as GL outer layer cloak and the 
GL inner layer cloak (3)- (5) can be coupled to construct 
new GLLH EM double layer cloak. 



III. GL EM CLOAK INVERSION 



B. FUNCTIONS IN GL EM CLOAK INVERSION 

We propose a GL EM cloak inversion for invisible cloak 
in this section. We request two invisible functions and 
without exceed light speed violation for the GL EM cloak: 
(I) the exterior EM field should not be scattering inter- 
fered by the cloak ; (II) the exterior EM field can not 
penetrate into the concealment of the cloak; and (HI) 
without exceed light speed violation. 



We proposal a GL EM inversion for cloak in this sec- 
tion. The 3D EM integral eqution and GL EM modeling 
in [2] and GL Metro Carlo inversion in [3] are base of the 
GL EM cloak inversion method. 



C. GL EM CLOAK EXTERIOR INVERSION 

When the source r^ and observer r are located outside 
of the cloak, we present the cloak invisible function I as 



the following inversion. 



GL EM CLOAK INVERSION 



Exterior EM field must not be interfered by scattering 
from the cloak. 

According to the invisible function I (6), and using the 
EM integral equation (8), we have 
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Using the EM integral equation (9), we have 



= 0, 



(11) 



The EM integral equation (11) and (12) for rg > Ro, 
and r > Ro, the EM integral equation (13) and (14) for 
rs > Ro, and r < Ri, and the "No exceed light speed" 
radial dependent constraint (15) are coupled to construct 
the GL EM cloak inversion for the EM invisible cloak. We 
use the GL EM Metro Carlo inversion method in [3] to 
solve the GLLH EM cloak inversion, we find the radial 
dependent GLLH EM cloak material to satisfy the GL 
EM cloak exterior inversion (11) and(12); the GL EM 
cloak inner inversion (13) and (14); and constraint (15) 
without exceed light violation. 
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D. GL EM CLOAK INNER INVERSION 



IV. NOVEL EM PROPAGATION IN GLLH EM 
= 0, CLOAK WITHOUT EXCEED LIGHT SPEED 

(12) In this section, we present a novel EM propagation in 

the GLLH EM cloak without exceed fight speed. The EM 
wave proagation pattern completely is new, and never be 
shown in other authors' papers. 



When the source r^ is located outside of the cloak and 
observer r is inside of the cloak , we present the cloak 
invisible function II in the following inversion. 

Exterior EM field does not penetrate into concealment. 

According to the invisible function II (7) and using the 
EM integral equation (8), we have 



Eb{r,rs,t) 
Hb{r,rs,t) 
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E{r',r,,t) 
H{r',r,,t) 

(13) 
By using the EM integral equation (9), we have 



A. THE SIMULATION MODEL OF THE GLLH 
EM CLOAK 

The simulation model: the 3D domain is 
[— 0.5m, 0.5to] X [— 0.5m, 0.5?7i] x [— 0.5TO,0.5m], the 
mesh number is 201 x 201 x 201, the mesh size is 0.005m. 
The electric current point source is defined as 
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0, 

(14) 



E. CONSTRAINT IN GLLH EM CLOAK 
INVERSION 

We present the WITHOUT EXGEED LIGHT SPEED 
VIOLATION to be the constrainst of the GLLH EM 
cloak inversion, for the radial dependent e(r) and p-{r), 



er{r)tie{r) > 1, 
ee{r)iir{r) > 1, 
eeir) = e^{r) = fig{r) = ^^(r). 



(15) 



S{r - r,)6{t)e, 



(16) 
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FIG. 5: (color online) At time step 40di, front of 

Electric wave, E^x inside of the GLLH EM cloak R\ <t < R2 
propagates no faster than light speed. The red S in right is 
source. 



where the r^ denotes the location of the point source, 
the unit vector e is the polarization direction, the time 
step dt = 0.3333 x 10^^° second, the frequency band is 
from 0.05 GHz to 15 GHz, the largest frequency / = 
15GHz, the shortest wave length is 0.02m. The GLLH 
EM cloak is consist of the spherical annular Ri < r < R2 
with center in original point and is fill in by the GLLH 
EM cloak materials , where i?i = 0.2491 meter, R2 — 
0.47 meter. The cloak is divided into 90 x 180 x 90 cells. 
The spherical coordinate is used in the sphere r < R2, 
the Cartesian rectangular coordinate is used in outside 
the sphere r > R2 to mesh the domain. 



B. NOVEL ELECTRIC WAVE PROPAGATION 

The electric intensity wave is excited by the point 
source S, which is denoted by red S in the figures in 
this paper. In Figures 1-2 and Figures 5-18, the source 
by red S is located in free space, in the right side outside 
of the whole GLLH EM cloak, at (0.83to, 0.0, 0.0). In 
Figure 19, the point source by red S is located in the left 
outer side of the cloak, at (— 0.83r7i, 0.0, 0.0). In Figure 
20, the point source by red S is located in the bottom 
outer side of the cloak, at (0.0, — 0.83to, 0.0). 

In Figure 5, at time step AOdt, electric wave, E^x inside 
of the GLLH EM cloak Ri < r < R2 propagates no faster 
than light speed. Its front inside of GLLH EM cloak is 
behind of the electric wave front in free space 

The Figure 6 shows that at time step 50dt, 
electric wave, E^x inside of the GLEM cloak Ri < r < 
R2 propagates no faster than light speed. 



In Figure 7, at time step QQdt, eElectric wave, Exx 
inside of the GLLH EM cloak Ri < r < R2 propagates 
slower than light speed. Its front inside of GLEM cloak 
is behind of the electric wave front in free space 

The Figure 8 shows that at time step 70dt, 
electric wave, Exx inside of the GLLH EM cloak Ri < 
r < R2 propagates slower than light speed. 

In Figure 9, at time step 80dt, electric wave, Exx inside 
of the GLEM cloak Ri < r < R2 propagates slower than 
light speed. Its front inside of GLLH EM cloak is behind 
of the electric wave front in free space 

The Figure 10 shows that at time step 90dt, 
electric wave, Exx inside of the GLEM cloak Ri < r < 
i?2 propagates slower than light speed. 

In Figure 11, at time step lOOdt, electric wave, Exx 
inside of the GLLH EM cloak Ri < r < R2 propagates 
slower than light speed. Its front inside of GLEM cloak 
is behind of the electric wave front in free space. 

The Figure 12 shows that at time step llQdt, 
electric wave, Exx inside of the GLLH EM cloak Ri < 
r < R2 propagates slower than light speed. 

In Figure 1 and Figure 13, at time step ll9dt, 
electric wave, Exx inside of the GLEM cloak Ri < r < 
i?2 propagates slower than light speed. Its front inside of 
GLLH EM cloak is behind of the electric wave front in 
free space. 

The Figure 2 and Figure 14 shows that at time step 
120dt, electric wave, Exx inside of the GLEM cloak 
Ri < r < R2 propagates slower than light speed. The up- 
side and downside parts of curved electric wave front are 
intersected at a branching point. These branching points 
form a 2D subsurface which depends on the source loca- 
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FIG. 6: (color online) At time step 50dt, front of 

Electric wave, Exx inside of the GLLH EM cloak R\ <r < R2 
propagates no faster than light speed, the red S in right is 
source. 



FIG. 7: (color online) At time step 60di, front of 

Electric wave, Exx inside of the GLLH EM cloak R\ <r < R2 
propagates slower than light speed, the red S in right is 
source. 
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FIG. 8: (color online) At time step 70dt, front of 

Electric wave, Exx inside of the GLLH EM cloak R\ <r < R2 
propagates slower than light speed 



FIG. 10: (color online) At time step 90dt, front of 
Electric wave, E^x inside of the GLLH EM cloak Ri < r < R2 
propagates slower than light speed 
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FIG. 9: (color online) At time step 80dt, front of 

Electric wave, Exx inside of the GLLH EM cloak R\ <r < R2 
propagates slower than light speed 



FIG. 11: (color online) At time step lOOdt, front of 
Electric wave, Exx inside of the GLLH EM cloak Ri < r < R2 
propagates slower than light speed 



tion. In this Figure and Figure 2, the red S denotes the 
source which is located in the right of the cloak. The front 
branching point is located in left of the concealment. In 
Figure 15, at time step 121dt, electric wave, E^x inside 
of the GLLH EM cloak Ri < r < R2 propagates slower 
than light speed. At time step 121dt, the wave front in- 
side of GLEM cloak is split into two wave fronts at the 
branching point. One wave front continously propagates 
outgoing and slower than the light speed. Other closed 



wave front propagates attracting into the inner bound- 
ary r — Ri oi the cloak and much slower than the light 
speed. The speed of the attracting wave front is going to 
zero. Its amptulide is rapidly decay to zero. 

In Figure 16, at time step 122dt, electric wave, E^x 
inside of the GLLH EM cloak Ri < r < R2 propagates 
slower than light speed. Its outgoing wave front inside of 
GLLH EM cloak is behind of the electric wave front in 
free space. Other attracting wave front propagates much 
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FIG. 12: (color online) At time step llOdt, front of 
Electric wave, Exx inside of the GLLH EM cloak Ri < r < R2 
propagates slower than light speed 
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FIG. 13: (color online) At time step llSdt, front of 
Electric wave, Exx inside of the GLWF cloak Ri < r < R2 
propagates slower than light speed 



slow than the light speed. Its speed is going to zero. Its 
amptulide is rapidly decay to zero. 

In Figure 17, at time step 130di, electric wave, E^x 
inside of the GLEM cloak Ri < r < R2 propagates slower 
than light speed. Its outgoing wave front inside of GLLH 
EM cloak is behind of the electric wave front in free space. 
Other attracting wave front propagates much slow than 
the light speed. Its speed is going to zero. 

In Figure 18, at time step 133dt, the outgoing front 
of electric wave, Exx propagates out of the GLEM cloak 



FIG. 14: (color online) At time step 119dt, front of 
Electric wave, Exx inside of the GLLH EM cloak R\ <r < R2 
propagates slower than light speed. The curved front inter- 
sects at branching point in left which depends on the source. 
The red S in right is source. 
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FIG. 15: (color online) At time step 120dt, front of 
Electric wave, Exx inside of the GLLH EM cloak R\ <r < R2 
propagates slower than light speed. The curved front is split 
to two branching fronts 



i?i < r < i?2 and recovers to wave in free space and is 
not scattering interfered by the cloak Other attracting 
wave front propagates much slow than the light speed. 
It propagates going to the inner boundary r — i?i of the 
cloak. Its speed is going to zero. 

In the Figure 1 and Figure 2, and Figure 5-18, the 
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FIG. 16: (color online) At time step 122dt, front of 
Electric wave, Exx inside of the propagates slower than light 
speed 



FIG. 18: (color online) At time step 133dt, front of 
Electric wave, Exx inside of the GLLH EM cloak Ri < r < R2 
propagates slower than light speed 
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FIG. 17: (color online) At time step 130dt, front of 
Electric wave, Exx inside of the GLLH EM cloak Ri < r < R2 
propagates slower than light speed 



FIG. 19: (color online) At time step 119dt, front of 
Electric wave, Exx inside of the GLLH EM cloak R\ <r < R2 
propagates slower than light speed. The red S in left is source. 



source is denoted by the red S which is located in the 
right of the cloak. 

Figure 19 shows that at time step 120dt, electric wave, 
Exx inside of the GLLH EM cloak Ri < r < R2 propa- 
gates slower than light speed. The upside and downside 
parts of curved electric wave fronts are intersected at a 
branching point in right. These branching points form 
a 2D subsurface which depends on the source location. 
In this Figure, the red S denotes the source which is lo- 



cated in the left of the cloak. The front branching point 
is located in right of the concealment. 

In Figure 20, at time step 120di, electric wave, E^x 
inside of the GLEM cloak Ri < r < R2 propagates slower 
than light speed. The left and right parts of the curved 
electric wave fronts are intersected at a branching point 
in topside. These branching points form a 2D subsurface 
which depends on the source location. In this Figure, the 
red S denotes the source which is located in the bottom of 
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FIG. 20: (color online) At time step 119dt, front of 
Electric wave, Exx inside of the GLLH EM cloak R\ <r < R2 
propagates slower than light speed. The red S in bottom is 
source. 



the cloak r — Ro, e and fi in (2)- (5) satisfy 

Sr == £0 = £0 = 1, re Ro, 
Mr = /^e = M0 = 1) re Ro- 



(17) 



This is necessary condition for the invisibihty of the 

cloak. From the GLLH EM cloak exterior inversion 

(11) and (12), we can prove that the condition (17) is 
necessary. 

Statement 2 Any exterior EM wave field can not 
penetrate into the concealment. 

From the GLLH EM inner cloak inversion (13) and 
(14), and EM integral equation (8) and (9), the state- 
ment 2 is proved. Also, because on the inner boundary 
r = i?i,the refractive index yf£rlJ^9 or yeg/Tr in (4)-(5) is 
infinite, from the convergence of integral (11)-(14), the 
EM wave field must be vanished on the inner boundary 
r = Ri. We obtain the condition 



Eg = E^ = Q, reR^ 
Hg = H^^O, re R, 



(18) 



Inversely, the conditions (17) and (18) are sufficient con- 
dition for the invisibility of the GLLH EM cloak in . 



the cloak. The front branching point is located in topside 
of the concealment. 



V. INVISIBLE PROPERTIES OF GLLH EM 
CLOAK 

By using the GL EM Metro Carlo inversion method 
in [3], we solve the GL EM cloak inversion (11), (12), 
(13), (14) and (15), and obtain the GLLH EM invisible 
cloak . The invisible properties of the GLLH EM cloak 
are verified in this section. 



A. 1. The Invisibility of The GLEM Cloak 

Statement 1 The exterior EM field is not scattering 
interfered by the cloak 

Because our GLLH EM exterior cloak inversion, by 
substituting the GLLH EM cloak in (2)-(5) and EM wave 
field by GL modeling [2] into (11) and (12), such that the 
integral equation (11) and (12) were held, we have iden- 
tity (6) for any exterior source and frequency, therefore, 
any exterior EM wave field must not be scattering in- 
terfered by the cloak. By the 3D GL EM modeling [2] 
simulations, electric wave propagations in the figures 5-18 
verify that there is no any reflction and no any scatter- 
ing from the cloak to interfer the exterior EM field in free 
space. 

By calculations, when Ri — 0.2491 meter and Ro = 
0.47 meter, we can prove that on the outer boundary of 



B. Without Exceed Light Speed Violation 

Statement 3 The refractive index in the GLLH EM 

cloak configration in this paper is large than one, the 
group speed in the cloak is finite and less than one. 

From (4) , the Sr and fir are finite positive and monotone 
increase function on variable r in the Ri < r < Ro- 
moreover, their minimum on the inner boundary of the 
cloak are 



iR,) = fir{.Ri) 



2Rg{Ro — Ri) 

i?2(l+4i?o + 2i?2) 



> 0, (19) 



their maximum value on the outer boundary Ro = 0.47 
of the cloak are one. 

From (5), the eg and fig are monotone decrease func- 
tion on variable r in the Ri < r < Ro, moreover, their 
maximum value on the inner boundary of the cloak are 
infinite, their minimum value on the outer boundary, 
Ro = 0.47 meter of the cloak are one 

The refractive index of the GLLH EM cloak in (2)-(5) 
is. 



nir) = ^/erfJ.e = VmJ^ > 1- 



(20) 



The group speed is finite and Vg < I. Therefore, in the 
annular configuration Ri < r < Ro, with Ri = 0.24 and 
Ro = 0.47 The GLLH EM invisible cloak material in (2)- 
(5) has no violation of exceeding light speed. The novel 
electric intensity wave propagations, which are shown in 
Figure 5 -18, in the GLLH EM cloak in (2)-(4) verify that 
the GLLH EM cloak is invisible cloak without exceed 
light speed violation. 
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C. Wave Front Branching 



The EM wave propagation in the GLLH EM cloak is 
completely novel that are displaying in Figure 1 and Fig- 
ure 2 and in Figures from 10 to 19. In particular, in 
Figure 14, EM wave front is curved around the conceal- 
ment and its upside and downside two parts intersect at 
a moveable branching point, which form a 2D subsurface. 
The Figure 15 shows that the two wave fronts are sep- 
arately propagation from the branching point. The first 
wave front is forward propagating outgoing. The second 
wave front is attracting propagating and approaching to 
the inner boundary r = Ri. The speed of two wave front 
propagation are shower then the light speed. The sec- 
ond wave front attracting propagating is very shower. 
Its amplitude is rapidly decay to zero. Its propagation 
speed is going to zero. The front branching points form 
a moveable 2D subsurface which depend on the source 
location. The front branching point and source point are 
separately located on two sides of the concealment room. 
In Figure 14 and Figure 15, the source is located in the 
right side of the cloak which is denoted by red S, the front 
branching point is located at the left side of the conceal- 
ment. In the Figure 19, the source, which is denoted by 
red S, is located in left of the cloak, the front branch- 
ing point is located in the right side of the concealment. 
In the Figure 20, the source, which is denoted by red 
S, is located in bottom of the cloak, the front branching 
point is located in the top side of the concealment. These 
novel EM wave propagations in our GLLH EM cloak and 
wave front branching phnomena only display in our pa- 
per. They never appear in other cloaks and never appear 
in other materials. The GLLH EM cloak , novel wave 
propagations and wave front branching are patent and 
copyright by authors in GL Geophysical Laboratory. 
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FIG. 21: (color online) At time step 84di, front of 
Electric wave, E^x inside of the GLLH EM cloak R\ <r < R2 
propagates slower than light speed. The red S in right is 
source. 
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VI. GLLH EM DOUBLE LAYER CLOAK 



We propose new GLLH EM double layer cloak in this 
section. 

Statement 4 When the concealment room of the 
single layer cloak is fill in the basic EM materials Eq and 
/iQ or normal material with refractive index n > 1, and 
the cloak materials are finite, then there is no EM field 
can be excited by source inside of the concealment. 



FIG. 22: (color online) At time step 84di, front of 
Electric wave, E^x inside of the PS cloak Ri < r < R2 
propagation is split two phase around the concealment, the 
front phase exceed light speed, behind phase slower than light 
speed. The red S in righ is source. 



We have proved the statement in [10] and [15] . The 
statement 4 theoreticaly confirms that our double layer 
cloth phenomena in 2001 in [1] is a credible physi- 
cal discovering. We propose GL double layer cloak in 
[9][11],[12],[13], and [16] and new GLLH EM double cloak 
here without exceed light speed violation. The inner layer 
cloak of The GLLH EM double layer is same as in (1) of 
[11] or (1) of [9], [12], [13], and [16], The GLLH EM cloak 



12 



GLW EM CLjOM< for BROAD FSEQUEWCt HANO 
WTTHOUT E;(CEED LIGHT SPEED VFOLATiOW 

Msds Ev L^ji Xie and J:3nhua LI GXIE In GLGED 
Copyrijfit jnd Patent ira R sservefi By GLGEQ 




linn 



-D.:5 i3 61i' i).5 

Harizontal Distance X ^M}, KadB By GL GEO 



FIG. 23: (color online) At time step 9Qdt, front of 
Electric wave, Exx inside of the GLLH EM cloak Ri < r < R2 
propagates slower than light speed 
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FIG. 24: (color online) At time step 90di, front of 
Electric wave, Exx inside of the PS cloak Ri < r < R2 prop- 
agation is much more exceed light speed. The red S in righ is 
source. 



is used as GL outer layer cloak of new GLLH EM double 
layer cloak. In the GL inner cloak, the EM wave can 
be excited by source inside of the concealment, when the 
EM wave propagates enter the GL inner layer cloak, it 
is rapidly decay, very slower than light speed and for- 
ward going to zero at outer boundary of the inner cloak. 
There is no any reflection back to the concealment from 
the GL inner layer cloak. The EM environment of its 
concealment is normal. The GL outer layer cloak has 



full invisible functions for broad frequency band, has re- 
fractive index large than one, and has no exceed light 
speed violation. The reciprocal law is satisfied. There- 
fore, our GLLH EM double layer cloak has advantages to 
overcome the three physical violations in Ps cloak or in 
other cloak. The Ps cloak is single cloak with finite di- 
electric and permeability materials. By the statements, 
there is no EM wave can be excited by source inside con- 
cealment of the Ps cloak. So Ps cloak's third physical 
violation caused that its concealment is dark room. 



VII. EM WAVE FIELD DECAY TO ZERO 
INVERSE RADIAL 

Statement 5, When source is located outside of the 
GLLH invisible cloak in (2)- (5), the observer r is located 
inside of the cloak and is going to the inner boundary 
r = Ri, the EM wave field delay to zero inverse radial. 
We have 



Er{r,e, 
Hr{r,9, 



Ee{r,e,i 



-C/0 



Q,when r ^ Rt, , . 
- 0, when r — Ri, ^ ' 



where r = Ri is the inner boundary of the GLLH cloak, 
on the inner boundary r = Ri, the radial electric wave 
field Er{r,9,(j)) = and the radial magnetic wave field 
Hr{r, 9,(1)) = are distictive property of the GLLH invis- 
ible cloak. On the inner boundary of other transformed 
cloak, the radial wave EM field are not zero. 



VIII. HISTORY AND DISCUSSION 

We discovered the GILD double layer cloth phenom- 
ena [1] to prevent detection from the exterior EM wave 
in 2001. The double layer cloth was appeared in residual 
magnetic field Hy in the GILD EM modeling and inver- 
sion [1] . Many times repeat GILD simulations show that 
the GILD double layer cloth is a mathematical physical 
phenomenon. We develop complete new GL EM model- 
ing [2] and GL EM inversion [3] [courant] in 2003 to inves- 
tigate the strange phenomena. We did find mirage [5] and 
double layer cloaks [9-13]. Our GLLH EM cloaks in this 
paper and in previous papers are based on the GL EM 
modeling and inversion and number theory. We deeply 
investigate the relationship between the field and mate- 
rial interactive scattering and non scattering and propose 
the GLLH EM cloak inversion. Our GL methods do not 
need artificial boundary condition to truncate the infinite 
domain. The GL method does not need to solve big ma- 
trix equation. The GL method combines the analytical 
method and numerical method. It can perform analyti- 
cal and numerical or mixed field simulation and material 
generation. Front branching is complete novel propaga- 
tion. Attracting front has benefit for detect the exterior 
field information in invisible cloaked concealment. 

The GLLH EM cloak materials (4)-(5) is similar with 
the form N"" / log'^ N in number theory [20] [21] [22] and 
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/i" log^(l/ft.) in superconvergent estimation [5][23]. That 
is consistent with our double layer cloth phenomena dis- 
covered in the residual field in [1]. By paper [20] and 
[21], GOLDBACH's conjecture proof will have a big step 
forwarding. 

The electromagnetic physics is governed by Maxwell 
equation. The equation is main tool for physics that 
relative both of geometry and number theory. Our GL 
EM full wave field modeling and inversion are based 
on the Maxwell equation The form N^/log^N and 
h°'\og' {1/h) are important mathematical form in the 
number theory and numerical estimation A good luck 
thing is that we chose these form and their appreciate 
combination and GL EM inversion to discover the invis- 
ible cloak material in without exceed light speed. Other 
cloaks based on the coordinate geometry transformation 
have exceed light speed and infinite speed physical viola- 
tions. 

To compare with Ps cloak which is based on the coor- 
dinate transformation, we simulate EM full wave propa- 
gation through the Ps cloak material using our GL EM 
modeling with same source, frequency band and domain 
geometry as our GLLH EM cloak simulation. In Fig- 
ure 21, at 84 time step, the EM wave propagates in the 
GLLH EM cloak and slower than the light speed. Its 
front inside the GLLH EM cloak is behind of the front 
in the free space. In Figure 22, at 84 time steps, the 
EM wave propagates through Ps cloak, its front is split 
to two phases, its front phase exceed the light speed. At 
90 time steps, the Figure 24 shows that the speed of EM 
wave inside Ps cloak much more exceed the light speed. 
At 90 time step, the EM propagation in GLLH EM cloak 
is presented in Figure 23, it is obvious that the EM wave 
front inside GLLH EM cloak is behind of its front in the 
free space which shows that the speed of EM wave inside 
GLLH EM cloak is not exceeding light speed. We propose 
GLLH EM cloak without exceed light speed in ellipsoid, 
2 polar geometry domain and other geometry domain. 
For overcoming the infinite speed in [14], Ulf Lconhardt 
et al. proposed a new cloak with finite speed based , on 
a Euclid and non Euclid joint transform[17]. The ULF 
cloak overcomes the "infinite speed" physical violation, 
even though its refractive index less than one in some 
subdomain, for example, n{(T,a',T) — n(0.757r, 7r,r) < 1. 
Our GLLH EM cloak without exceed light speed in the 
2 poplars geometry domain and in the ellipsoid domain 
will be published. 



IX. CONCLUSION 

The GL EM modeling and inversion theoretical analy- 
sis and many GLLH EM modeling simulations verified 



that the GLLH EM cloak has full invisible functions. 
The GLLH EM invisible cloaks have no infinite speed 
and exceed light speed two violations. The novel EM 
wave propagations have presented in Figure 1-2 and Fig- 
ure 5-18. The EM wave front inside of the GLLH EM 
cloak has always been behind of the wave front in the 
free space. The front inside of the GLLH EM cloak 
completely does curve around the concealment and form 
curved Grescent-like in Figure 1 and 13. The two front 
teeths of the Orescent front intersect at the branching 
point in Figure 14, 19 and 20. The branching points form 
2D subsurface. The wave front is split to two fronts at 
the front branching point. The outgoing front propagates 
forward outgoing and recovers to original wave front in 
free space. The attractive front propagates and to the 
shrinks to inner boundary of the cloak. Its amplitude 
and speed are rapidly decay to zero. Using GLLH EM 
cloak as outer layer, the GLLH EM double cloak has nor- 
mal environment concealment. In its concealment, the 
EM wave field excited by the internal source propagates 
enter to the inner layer cloak, its amplitude and speed 
rapidly decay to zero before the outer boundary of the 
inner layer cloak. Therefore, the GLLH EM double invis- 
ible cloak overcomes the three physical violations in other 
cloaks. GLLH EM double invisible cloak can be practica- 
ble. GLLH EM cloak software can generate various scale 
full invisible cloaks without exceed light speed, and can 
generate GLLH EM double cloak for various scale that 
have no physical violation and have wide applications. 
GLLH EM cloak and GLLH EM cloak modeling and in- 
version and GLLH EM cloak software are patented and 
copyrighted by authors in GL Geophysical Laboratory 

The GL modeling and inversion method is an effec- 
tive physical simulation method. It has double abilities 
of the theoretical analysis and numerical simulations to 
study the cloak metamaterials and wide material and 
Field scattering in physical sciences. 
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